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Abstract 

The sign-reversing fully gapped superconducting state, which is expected to be realized in oxyp- 
nictide superconductors, can be prominently affected by nonmagnetic impurities due to the inter- 
band scattering of Cooper pairs. We study this problem based on the isotropic two-band BCS 
model: In oxypnictide superconductors, the interband impurity scattering /' is not equal to the 
intraband one /. In the Born scattering regime, the reduction in is sizable and the impurity- 
induced density of states (DOS) is prominent if / ~ /', due to the interband scattering. Although 
impurity-induced DOS can yield a power-law temperature dependence in 1/Ti, a sizable suppres- 
sion in Tc is inevitably accompanied. In the unitary scattering regime, in contrast, impurity effect 
is very small for both Tc and DOS except at / = /'. By comparing theory and experiments, we 
expect that the degree of anisotropy in the s±-wave gap function strongly depends on compounds. 



I. INTRODUCTION 



Recently, the mechanism of superconductivity in high-Tc superconductors with FeAs lay- 
ers [1-5] has been attracting considerable attentions. The superconducting state is reahzed 
by introducing carrier into the parent compound, which shows the spin density wave (SDW) 
state at ~ 130K [6, 7]. In the SDW state, the ordered magnetic moment is ~ 0.3 /xb and 
the ordering vector is Q (tt, 0) [6, 8-10]. NMR studies had clearly shown that the singlet 
superconducting state is realized in iron oxypnictides [11-13]. A fully gapped superconduct- 
ing state has been determined by the penetration depth measurement [14], angle- resolved 
photoemission spectroscopy (ARPES) [15-18], specific heat measurement [19], and so on. 

In the first-principle band calculations [20-22], the Fermi surfaces in iron oxypnictides are 
composed of two hole-hke Fermi pockets around the F = (0, 0) point and two electron-like 
Fermi pockets around M= (tt, 0), (0, tt) points. The nesting between the hole and electron 
pockets is expected to give rise to the SDW state in undoped compounds. In doped com- 
pounds without SDW order, the antiferromagnetic (AF) fluctuations with Q fa (tt, 0) is 
expected to induce a fully gapped s-wave state with sign reversal, which is called the s±- 
wave state [23-30]. Moreover, near the SDW boundary, the Hall coefficient and Nernst signal 
show prominent anomalous behaviors [7, 12, 31], which are similar to those observed in high- 
Tc cuprates and in CeMIns (M=Co,Rh,Ir) [32]. Theoretically, these anomalous transport 
phenomena indicate the existence of strong AF fluctuations [33]. At the same time, huge 
residual resistivity far beyond the s-wave unitary scattering is expected to appear near the 
SDW boundary theoretically [34] . 

To investigate the pairing symmetry of superconductivity, impurity effects on the su- 
perconducting state offer us decisive informations. In iron oxypnictide superconductors, 
impurity eff'ect on due to Co, Ni, or Zn substitution for Fe sites is very small or absent 
[11, 12, 35-37]. This result clearly rules out the possibility of line-node superconductivity. 
One may also expect that the s-wave state with sign reversal is also eliminated, since the 
Cooper pair is destroyed by the interband scattering induced by impurities. However, we 
have recently shown that T^. is almost unchanged by strong (unitary) impurities, since the 
interband impurity scattering potential /' is different from the intraband one / [38]. The 
reason for this unexpected result is that the effective interband scattering is renormalized to 
zero in the unitary limit except at / = /' in the T-matrix approximation. Therefore, the ex- 
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perimental absence of impurity effect on Tc in iron oxypnictides is well understood in terms 
of the s±-wave state. On the other hand, Tc will be prominently reduced by short-range 
weak (Born) impurities [38]. 

Recently, several authors had revealed that in-gap density of states (DOS) is induced by 
impurities in the s-|--wave state using the Born approximation for general value of a; = \I' /I\ 
[39], or using the T-matrix approximation only for a; = 1 [40, 41]. They also demonstrated 
that the relation 1/71 ocT^ under Tc, which had been reported by several groups [13, 42-44], 
can be reproduced by the impurity-induced DOS. However, the assumed impurity parameters 
(nimp, I and /') also yields a sizable suppression in Tc according to the analysis in Ref. [38]. 
Furthermore, impurity- induced DOS should be sensitive to the vahic of x in the unitary 
scattering regime, as suggested in ref. [38]. Therefore, we have to study the impurity effects 
on the DOS and Tc for general x, and compare their relationships in detail. 

In this paper, we investigate the impurity-induced DOS and Tc in the s±-wave state using 
the T-matrix approximation for general x. We stress that x is not unity in iron oxypnictides 
since hole and electron pockets are not composed of the same d-orbitals. In the Born or 
intermediate scattering regime, a sizable impurity- induced DOS appears for x > 0.7, and 
therefore 1/Ti may deviate from a simple exponential behavior. Although impurity-induced 
DOS can yield a power-law temperature dependence in 1/Ti [39-41], we find that a sizable 
suppression in Tc is inevitably accompanied. The anisotropy in the s±-wave superconducting 
gap, which had been predicted theoretically [23, 25] , might be responsible for the power-law 
temperature dependence of 1/Ti under Tc as discussed in ref. [45]. In contrast, unitary 
impurities affect both the superconducting DOS and Tc only slightly, except at 7 = 

II. T-MATRIX APPROXIMATION IN THE TWO-BAND BCS MODEL 

As studied in refs. [29, 38, 39, 41], the s±-wave state is realized in the two-band BCS 
model if we introduce the interband repulsive interaction, which represents the AF fluctua- 
tions due to the interband nesting in iron oxypnictides. In the present paper, we study the 
impurity effect using the T-matrix approximation for general /'//. In the presence of mass 

enhancement due to many-body effect, m* /itlq > 1, both the superconducting gap and the 
impurity effect (or impurity concentration nimp) are renormalized by the factor (m*/mo)~^. 
In the present analysis, we neglect the mass-enhancement for simplicity. 
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In the Nambu representation, the two-band BCS model is given by [46, 47] 



(1) 



where cj^. 



(Ck},c^{,c_kiQ;,c^ki)> and 
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In eq. (2), e^, are the band dispersions measured from the Fermi level. Since we consider 
the isotropic s± superconducting state, only the DOSs for both bands at the Fermi level 
(A^Q,, Np) are taken into consideration in the present BCS study. Aq, in eq. (2) are the 
superconducting gap. When only the inter-band repulsive interaction {gap — g^a > 0) is 
taken into consideration, the gap equation without impurities is given as [38, 40, 41], 



MP) 



(3) 



where e„ = TrT{2n+l) is the fermion Matsubara frequency, and cUc is the cutoff energy. A^/3(q) 
is the DOS for /3(Q;)-band at the Fermi energy in the normal state per spin. //3(Q)(e) is the 
local anomalous Green function for j3{a) band, which will be given later. Since //3(a) oc Af3{a), 
the s±-state A„ = -Ap is realized for gais > [38, 40, 41]. Moreover, \Aa/Ap\ {N/s/Nay^'^ 
since fa/ ff3 ~ Aa/Ap. 

The Nambu matrix representation for the impurity potential is given as 

/ / J' \ 
r I 
0-7 -/' 

\ -r -I J 

We can assume that 1,1' > without losing generality. In the presence of impurities, the 
Green function in the Nambu representation is given by [47] 



I = 



(4) 



(5) 



where lo = cu + i5 {5 — +0), and S(a;) is the self-energy due to impurities. 



Hereafter, we derive S(a;) in the T-matrix approximation, which gives the exact result 
for riimp <S 1 for any strength of /, The T-matrix in the Nambu representation is given 

by 



T(u;)^(l-I-g(u;)y'l, 
where g{u;) = jj: ^j^, Gk(a') is the local Green function, which is given by [47] 

f g„{u) UCu) \ 
g^{u) fp{u) 

fa{Cb) ga{to) 

V fp{Cj) gf,{Cj).J 

In the above expression, g^ and /j (i = a, /5) are given by 

ujZAu) 



(6) 



(7) 



gi{ix!) = -nNi 



(8) 

(9) 
(10) 



where and E" (i = a, /?) are the normal and anomalous self-energies, respectively. In the 
T-matrix approximation, the self-energies are given by using eq. (6) as 



■?T'imp733(a;), E^(-a;) = -nimpT44(d;), 



(11) 
(12) 

(13) 



In the fully self-consistent T-matrix approximation, we have to solve eqs. (3) and (7)- 
(13) self-consistently. In this paper, however, we solve only eqs. (7)-(13) self-consistently, 
by neglecting the impurity effect on and Ag in eq. (3). This approximation is justified 
when the reduction in Aq,(^) due to impurity pair-breaking is small. 



III. NUMERICAL RESULTS 



Here, we discuss the impurity effect on the DOS and in the s±-wave superconducting 
state, based on the numerical results given by the T-matrix approximation. 
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FIG. 1: — ATc/riimp as a function of / given by the T-matrix approximation in the case of (a) 
Nq, = Np = 1 and (b) A^^^ = l^Np = 0.5. In both figures, the unit of energy is 1/Na, which 
corresponds to 18000 K for (a) and 14000 K for (b), since Na+Nfj = 1.31 eV~^ in iron oxypnictides. 
-ATc/nimp for (x = 1, / = oo) is l/SiV^ ~ 2300 K for (a), and l/5.84iV„ ~ 2400 K for (b). The 
superconducting DOS for parameters denoted by filled diamonds are shown in Figs. 3 and 5. 

A. Impurity effect on Tc 

As derived in Ref . [38] , the expression for the reduction in per impurity concentration 
based on the two-band BCS model is given as 

AT, _ 7r2[3(iV„ + iV^)-2v/AgV^]/'2 



For riimp -C 1, the transition temperature is given by = T° — (— ATc/riimp) ■ ^imp, where 
T° is the transition temperature without impurities. In eq. (14), A = 1 + it'^P{N^ + N^) + 
2NaNfsTT^r^ + N^N^TT\P-r^f. In the case of x = /'// = 1, the right hand side of eq. (14) 
is [3{Na+Nf3)-2y/N^N0]/[8{No,+Np)'^]+O{I-'^) in the unitary regime. In the case of x 7^ 1, 
in contrast, eq. (14) is given by x^[3{Na + N/s) - 2^/NJ^]/[87T^N^N^{1 - x^)/^] + 0{r^). 
Therefore, eq. (14) approaches zero in the case of x 7^ 1 in the unitary regime. 

Figure 1 (a) shows — ATc/nimp given in eq. (14) in the case of Na = Np = 1. In iron 
oxypnictides, the total DOS per Fe atom {Na + N^) is 1.31 eV~^ per spin [20]. Then, 
1/N = 1 corresponds to 18000 K. When x = 1, -ATjn^^p approaches 1/8A^ ~ 2300 K in 
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the unitary regime {IN ^ 1). Therefore, the superconductivity in iron oxypnictides will 
vanish only at riimp ~ 8N ■ T° = 0.01 ~ 0.02 [1 ~ 2 %]. When a; 7^ 1, in high contrast, 
— ATc/nimp decreases and approaches zero as I increases in the unitary regime, since the 
effective interband scattering is renormalized as Jg^ ~ /' • (IN)^"^ <^ I' [38]. 

According to the first principle calculations, N^/Na > 0.7 in iron oxypnictides [48]. Here, 
we study the case of Np/Na = 0.5 in order to clarify the the impurity effect on Tc for the the 
particle-hole asymmetric case; Np/Na 7^ 1. Figure 1 (b) shows — ATc/riimp for Na = 1 and 
Nf3 = 0.5. According to eq. (14), -ATju^^p for x = 1 and / = 00 is 1/5MN^ ~ 2400 K, 
by taking account of the relation Na + Nf^ = 1.31 eV~^ in iron oxypnictides. By comparing 
with the results for N^ = Njj = 1 in Fig. 1 (a), we find that — ATc/riimp is insensitive to the 
value of Np/Na, under the condition that Na + Np =constant. 

I' I' I' I' r 
— \ ( I y X + = ^0 

~ op OpO^ 1 + (7I/W)- 

apaapapa l + C'^^ ^ ) 

FIG. 2: Intraband and interband T-matrices in the normal state, Taa^ and T^J*^. 

Previously, impurity effect on Tc in two-band BCS models had been studied by many 
authors in various contexts [49-54], and it was found that Tc is unchanged in the unitary 
limit [50, 51]. However, eq. (14) for s±-state had not been derived. Here, we present a clear 
explanation why the interband scattering (pair breaking) is absent in the unitary regime, 
which had not been discussed previously. Figure 2 shows the intraband and interband T- 
matrices in the normal state, T^^" and T^J°, in the case of / = and Na = Np = N. 
Apparently, T^J° approaches zero for /' 00. Next, we consider Tajj for general {1,1'). If 
we construct Tap of {T^^^, T^J^, I), it contains at least one T^J°. For this reason, interband 
T-matrix is expected to approach zero in the unitary regime. This expectation is correct 
unless a; = 1, as shown in Ref. [38]. 



7 



(a) 

1=0.25 



(b) 

1=0.5 





FIG. 3: Obtained DOS in the superconducting state for Na = Np = 1 ('^ 0.66 eV ) and Aq, = 
-A^ = 0.005 (~ 100 K), in the case of (a) / = 0.25, (b) / = 0.5, (c) 1 = 2, and (d) / = 8. Impurity 
concentration riimp is 0.008. The insets in (b)-(d) present the DOS for n\yap = 0.001. Note that 
N{-uj) = N{uj). 

B. Impurity effect on the DOS 

In the s±-wave superconducting state, impurity interband scattering not only reduces T^, 
but also induces the in-gap state in the superconducting DOS [40, 49, 51, 55]. The DOS is 
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given by the imaginary part of the local Green function, which is expressed in eq. (8), as 
follows: 

N{uj) = --lm{gaiu;) + gp{u;)}. (15) 

TT 

If riimp <^ 1, the obtained DOS will be reliable for any / and /' in the present T-matrix 
approximation. Figure 3 shows the DOS in the superconducting state in the case of Na — 
N/s ^ 1 and A« = -Ap = 0.005 for riimp = 0.008. \Aa,p\ = 0.005 corresponds to 90 K. 
Experimentally, in Bao.6Ko.4Fe2As2, |A| = 11 ~ 12 meV for a Fermi surface (hole-like) and 
for 7 and 6 Fermi surfaces (electron-like), and |A| = 5.8 meV for f3 Fermi surface (hole-like) 
[18]. As shown in Fig. 3 (a), the superconducting gap is almost filled by the impurity- 
induced DOS when I — I' — 0.25 (Born regime), which corresponds to ~ 5000 K. In this 
case, the nuclear relaxation ratio 1/Ti shows a power- law temperature dependence since the 
impurity- induced DOS is approximately linear in cu, like in line- node superconductors. 

Figure 4 shows 1/Ti below for 1^1'^ 0.25, where a = {-KNIf/il + {-KNIf) = 0.38. 
We put |A«(^)| = COOSv^l - T/Tc and % = 0.002. The method of calculation is explained 
in Refs. [40, 55]. For x = 1.0 and 0.9, 1/Ti is inside of T^- and T^-lines for > T > O.lTc, 
consistently with the analysis by Parker et al. for a — 0.4 [40]. In these cases, however, 
reduction in T^. due to impurities, which is given by njmp times — ATc/riimp in Fig. 1 (a), 
reaches 13 K. The estimated reduction in Tc would be underestimated since — ATc/riimp is an 
increase function of riimp for Tc < T'°/2 [41, 47]. In all cases we have studied (Fig. 3 (a)-(d)), 
power-law behavior in 1/Ti for T <^ due to the galpess superconducting state always 
accompanies a sizable suppression in Tc, — ATc > 10 K, for |A| = 90 K. When |A| = 40 
K, the gapless superconducting state can be realized when —ATc ~ 6 K. Thus, it will be 
difficult to ascribe the experimental relation 1/Ti oc T^ below Tc [42-44] in clean samples 
with high Tc to the impurity effect. 

In the intermediate (/ = 0.5) or unitary (/ > 2) regime, in Figs. 3 (b)-(d), a large 
impurity- induced DOS appears at the zero energy in the case oi x = 1, which is consistent 
with previous theoretical studies [40, 55]. In this case, however, —ATc ~ 13 K for riimp = 
0.008. If we put X < 0.9, in contrast, —ATc in the unitary regime is prominently reduced 
as shown in Fig. 1. At the same time, the impurity-induced DOS quickly moves to the gap 
edge and disappears, as demonstrated in Figs. 3 (c) and (d). The reason for these results 
is that the effective interband scattering is renormalized as 1'^^ ~ /' • {IN)~'^ <^ I' in the 
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FIG. 4: Obtained 1/Ti below Tc for / = 0.25 and x = 0.7 ~ 1.0. For x < 0.8, 1/Ti decreases much 
faster than at low temperatures because of the absence of the m-gap state near w = 0. Since 
the quasiparticle damping rate is 7 = nimp7rA^(I^ + 1''^)/ A [38], 7/A = 0.18 for x = 1. 

unitary regime [38]. 

Next, we study the case of N^, = l,Np = 0.5. We put = 0.005 and = -0.0071 
by considering the relationship |Aq,/A^| ~ {Nf^/NaY^"^ in the two-band BCS model with 
repulsive interband interaction, as explained in §11. Figure 5 (a), (b) show the DOS in the 
superconducting state for riimp = 0.008. When I = I' = 0.25, in Fig. 5 (a), the impurity- 
induced DOS is reduced by changing A'^ from 1 to 0.5, by comparing with Fig. 3 (a). As x 
decreases from unity, the impurity- induced DOS moves to the gap edge. In the intermediate 
(/ = 0.5) or unitary (/ > 2) regime, impurity-induced DOS covers the zero energy state in 
the case of x = 1 and njmp = 0.008, as shown in Figs. 5 (b)-(d). However, finite gap appears 
around the Fermi level for a; < 0.9. In the case of / = 8, in Fig. 5 (d), impurity- induced 
DOS is very large at x = 1, whereas it is strongly suppressed for x < 0.9. When the impurity 
concentration is very low (riimp ~ 0.001), in-gap state deviates from uj = even if x = 1 as 
shown in insets in Figs. 5 (b)-(d), since fa + fi3 given in eq. (7) is non-zero in the case of 
Na^N^ [41]. 
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FIG. 5: Obtained DOS in the superconducting state for A^^^^ = 1,N^ = 0.5 and Aq, = 0.005, = 
—0.0071, in the case of (a) / = 0.25, (b) / = 0.5, (c) 1 = 2, and (d) 1 = 8. Impurity concentration 
nimp is 0.008. The insets show that impurity-induced DOS is always located at a finite energy for 
nimp = 0.001. 

IV. DISCUSSION 

In the present paper, we studied the impurity effects on the s-t-wave superconducting 
state, which is expected to be realized in iron oxypnictide superconductors. There, nonmag- 
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netic impurities can induce both the in-gap bound state and the reduction in T^. Based on 
the two-band BCS model, we have found that the zero-energy in-gap state emerge under 
the conditions that (i) x = \r/I\ = 1 and (ii) \I\Na,\I\Np » 1. Deviating from these 
conditions, in-gap state shifts to a finite energy, and disappears eventually. 

Here, we discuss the case of unitary scattering: In iron oxypnictide superconductors, Fe 
substitution by other elements (such as Co, Ni, and Zn) will cause the unitary scattering 
potential. In this case, the impurity potential is diagonal with respect to the d-orbital 
[38] . The impurity potential has off-diagonal elements in the band-diagonal representation. 
As discussed in Ref. [38], x ~ {Y.dC>d,a{^)0d,f3(k'))lla,k'eP^ where Od,«(k) = ((i;k|a;k) 
represents the transformation matrix between the orbital representation (orbital d) and the 
band-diagonal representation (band a). In iron oxypnictide superconductors, x ~ 0.5 since 
the hole-pockets are composed of dxz, dyz orbitals of Fe in the two-Fe unit cell, whereas half 
of the electron-pockets are composed of d^'^-y^ orbitals [38] . Since the impurity effect is weak 
except when x — 1 in the unitary regime as shown in Figs. 3 (d) and 5 (d), Fe substitution 
by other elements will affect the superconducting DOS and Tc only shghtly. 

We also discuss the case of Born scattering due to "in-plane" weak random potential or 
disorder: As shown in Figs. 3 (a) and (b), the impurity effect is rather insensitive to x. 
Therefore, a broad impurity-induced in-gap state will emerge in the superconducting DOS, 
and a sizable reduction in Tf. occurs at the same time. Born impurity scattering will be also 
caused by "off-plane" impurities like the As substitution by other elements. In this case, 
the radius of impurity potential R for Fe sites will be about the unit cell length a. Then, 
the impurity scattering (k k') is restricted to |k — k'| < ~ 1/a. Since |k — k'| ^ vr/a 
in the interband scattering between electron-pockets and hole-pockets, /' should be much 
smaller than /. Therefore, the effect of off-plane impurities on the s±-wave state will be 
small since the relationship x <^ 1 is expected to be realized. 

In summary, in iron oxypnictide superconductors. Born or intermediate in-plane impu- 
rities cause prominent impurity effects since the s±-wave state is violated by the interband 
scattering. Only one percent Born impurities with x > 0.5 induce not only plenty of in- 
gap DOS, but also sizable reduction in T^. For this reason, relation l/T-i oc T"^ below Tc 
observed in clean LaFeAsOi^x-F^. [42, 44] and in LaFcAsOo.7 [43] samples, which would be 
almost absent from the impurity reduction in T^, cannot be explained by the present analysis 
based on the isotropic BCS model. Thus, anisotropy in the s±-wave superconducting gap 



12 



might be responsible for the relation 1/Ti oc [45]. Recently, rapid suppression in 1/Ti 
(oc T"; q; > 5) below Tc had been observed in a clean LaFeAsOo.gFo.i sample with Tc — 28 
K (^intrinsic Tc) [57]. This result is consistent with the penetration depth [14] and ARPES 
[15-18], and it is naturally explained by the present analysis. Theoretically, in fully gapped 
s-wave superconductor, the gap function becomes anisotropic due to magnetic fluctuations, 
in a way that the two superconducting gap minima are connected by the nesting vector [56] . 
In iron oxypnictides, the degree of anisotropy in the s±-wave gap function is rather sensitive 
to model parameters such as the nesting condition [23, 25] . The wide variety of behaviors 
in 1/Ti would reflect the large sample dependence of the gap anisotropy in iron oxypnictide 
sup er conduct ors . 
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